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We determine the N -particle stationary states of a staggered stochastic hop- 
ping model with reflective boundaries. It is shown that the stationary states 
are in fact so-called optimum ground states. Recursion relations in the parti- 
cle number for any l-point density correlation function will be derived. Fur- 
thermore, the connection between reflective boundaries and the occurrence of 
optimum ground states is examined. An explicit counterexample shows that 
reflective boundaries do not enforce the stationary state to be an optimum 
ground state. 
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1. INTRODUCTION 

Interacting many-particle systems encountered in nature are sometimes 
difficult to describe in terms of classical or quantum mechanics. A stochastic 
description is then more appropriate for systems which behave essentially 
randomly on a phenomenological level. Often this leads to problems far away 
from thermal equilibrium for which - in contrast to equilibrium systems - a 
general concept or framework is still missing. Therefore one has to rely on 
the investigation of special examples. Models which are exactly solvablellil 
become then even more important than in equilibrium physicJll. They allow 
to probe the nature of non-equilibrium physics with the problems related 
to the reliability of approximative methods. Also stochasjtiiL systems are 
often used for the modelling of interdisciplinary problemsHoQ'El where the 
interactions between the entities are not known exactly. 

The dynamics of stochastic systems is governed by the master equation 
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which can be interpreted as a Schrodinger equation in imaginary timJi (see 
Sec. |2]) . In the analytical treatment of non-equilibrium systems it is therefore 
possible to exploit this similarity by using successful techniques from equi- 
librium physics. Especially the so-called matrix product Ansatjlp fMPA). 
first developed for spin chains0, is important. It has been shownHIll that 
the stationary state of one-dimensional stochastic processes is generically 
of MPA-type. However, in most cases it is still very difficult to calculate 
expectation values explicitly. 

In this paper the exact stationary state of a one-dimensional stochastic 
process with staggered interaction and refiective boundaries is determined. 
This model can be regarded as the merger of the ?7g[S'f/(2)]-symmetric 
hopping model of Sandow and Schiitai^ and a ratchet model proposed by 
Kolomeisky and WidomllZI (for details, see t3). A similar stochastic hop- 
ping model with periodic boundary conditions, where also coagulation and 
decoagulation processes may occur, has been studied by Fujii using free 
fermion techniques. Our model can not be solved using free fermion tech- 
niques but requires a different approach. This approach relies heavily on 
the observat ion tha t the stationary state of our model is an optimum ground 
state (OGS) The occurrence of OGS in the context of matrix prod- 

uct states will be elucidated. It will be shown that OGS are not ubiquitous 
for systems with reflective boundaries. 



2. QUANTUM FORMALISM AND DEFINITION OF THE 

MODEL 

We consider stochastic systems on a one-dimensional lattice with L 
sites. Each cell j can either be empty or occupied by a particle. Its state is 
characterized by the variable aj which could e.g. be the occupation number. 
The model is then defined by its dynamics, namely by the transition rates 
w{a — > P) for a transition from a lattice configuration a = (ai,... 
to the configuration /3 = (/3i,... The probability P{a;t) to find the 

system at time t in the configuration a = (ai, ... , ql) is then given by the 
master equation 

^EMl = Y^w{P^ a)PiP;t) -^wia^ (3)P{a;t). (1) 

Introducing the state vector \P{t)) = X^q, -P(ai , • • • , ol; • • • 
the master equation can be written in more compact form 

l-^\P{t)) = -H\P{t)), (2) 
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where the matrix elements of the stochastic Hamiltonian H can be expressed 
through the transition rates w{(3 — > a). For processes with nearest-neighbour 
interaction and periodic boundary conditions, the stochastic Hamiltonian H 
can be written as a sum H = Yl^=i °f local Hamiltonians acting 

only on the sites j and j + 1- For other boundary conditions the interactions 
at the boundaries j = 1 and j = L are described by additional boundary 
operators hi and hi. 

The so-called quantum formalism described above emphasizes the simi- 
larity with a Schrodinger equation in imaginary time. The stochastic Hamil- 
tonian H, which in general is not hermitean, is often related to spin chaingBllfl. 
Therefore it is not surprising that many methods for the treatment of equilib- 
rium systems can be adoptedStJ, e.g. free fermion techniques, Bethe Ansatz, 
or the matrix product Ansatz. 

The stationary state |P) reached for t — > oo is determined by 

H \P) = 0, (3) 

i.e. it corresponds to a zero-energy ground state of the Hamiltonian H. 

In this paper we consider the following stochastic process: Particles with 
a hard core occupy the sites of a chain of length L. The odd sites of the 
chain belong to sublattice A and the even sites to sublattice B. The global 
stochastic Hamiltonian is an alternating sum of local two-site Hamiltonians 

H = Y,f^3 hfj+i + Yl ^^j+i ' (4) 

jeA jeB 

where h^j^i and h^j^i act non-trivially on sites j and j + 1 according to 

/O 0\ 

c -i 
-c i 

(5) 

The local two-site basis is ordered as follows: {TT)Ti)iT)ii} where j corre- 
sponds to an empty site and J, to an occupied one. Here we have already 
used the spin language to emphasize the similarities with spin systems. The 
only allowed processes are jumps of particles to unoccupied neighbouring 
sites where the hopping rates alternate along the bonds of the chain. 

The /ij € are arbitrary constants which may differ from site to site. 
They control the activities of the bonds and influence dynamical properties 
only. 
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3. CONSTRUCTION OF THE STATIONARY STATE 

Since no particles enter or leave the chain, the total particle number 
is a conserved quantity and the stationary states of a chain of length L can 
be classified according to the particle number N < L. 

Exact solutions of small systems with L < 4 indicate that the stationary 
state for any given particle number is unique and, more interestingly, that 
this stationary state is an OGS. An OGS is a special type of ground state 
of a Hamiltonian H = Y2j i ^ji ' which is at the same time ground state 
of all local Hamiltonians hji. In the present case this means that all local 
Hamiltonians hfj^i and h^j-^-i have to annihilate the stationary state \P) 
since the "ground state energy" of is according to (^). Explicitly this 
has the consequence that if \P) has a component 

|ai,... ,aj_i, i, t,/3j+2, • • • ,I3l) 
it must also have the component 

|ai,... ,aj_i, t, i,/3j+2, • • • 
with a relative weight given by 

if j odd, 

-2 -f . W 
c II J even. 

Stochastic processes on a one-dimensional lattice with a finite interac- 
tion range and boundary interacti ons h five a stationary state which can be 
written as a matrix product state 113113. This can be shown without using 
the explicit form of the boundary operators hi, hi. Therefore, reflective 
boundaries^ - corresponding to vanishing boundary operators hi = = 
- are allowed as well. 

A suitable (grand canonical) matrix product ansatz for the stationary 
state in the presence of staggering has the form 



= \V), (7) 



L-fold product 

where we have used (^) to write the state for even and odd L in a unified 
way and 



D J ' 

^In the language of quantum spin chains such boundaries are usually denoted open 
boundaries. 
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E, D, E and D are matrices acting on some auxiliary vector space j4, \V) € A 
and (W\ G A*. By direct computation one verifies that this vector is a 
ground state of H, if one can find matrices Xi, X2 and Xi, X2, which are 
the components of the vectors X and ^, respectively, such that the following 
relations are fulfilled 

h^{V (S)V) = X(g)V-V0X, 

{W\X = 0, (9) 



X 



\V) = 0. 

L 



The most simple way to satisfy this algebra is given by the choice Xi 
X2 = Xi = X2 = 0. Then, the only remaining equations are 



h^{V®V) = 0, 
h^{t>0V) = 



B,^-^. . (10) 



turning |-Pl)o into an optimum ground state! 

More explicitly, the problem of finding the ground state of H becomes 
the problem of finding a representation of the quadratic algebra 

ED = irDE, 

f (11) 
ED = -^DE. 

The algebra ( pT|) fixes the relative weight of two configurations in \Pl)o that 
differ only by the interchange of a particle and a hole at the sites j and 
J + 1. This relative weight is exactly the one that has been proposed above 
by means of the optimum ground state property given by (^). Thus, in case 
of refiective boundaries the appearance of optimum ground states can be 
connected with the fact that the algebra necessary for the construction of 
the ground state in form of a matrix product state takes its most simple 
form. Nevertheless, it should be noticed that refiective boundaries do not 
enforce optimum ground states, as will be shown in the caveat at the end 
(Sec. i). 

Up to this point, two ways for the construction of the stationary state 
of H have been presented. The first way gives a recipe for the construction 
of the relative probabilities of all states that are present in the stationary 
state. The second way is connected with the matrix product approach: The 
ground state problem is transformed into the purely algebraic problem of 



finding the representation of the quadratic algebra (11). 



n 



(13) 
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In the following, no use of these two possibilities will be made. Instead, 
a third way will be given that closely res;mbles the construction in@. The 
starting point will be the vacuum vector |vac) := I T) describing a 

completely empty chain. 

Out of this vacuum all A^-particle stationary states will be constructed, 
using powers of suitable creation^ operators. With 

for every j G {1, 2, . . . ,L} a "local" particle creation operator bj is given by 

i=l 

The global creation operator B~ is the sum of all local creation operators 

L 

B- = Y.bJ. (14) 
i=i 

In order to shorten the notation one furthermore defines 

/ 1 \ l-^-l / 1 \ l-^J 1 - (ac)'' 

{k}a:c ■■= [-^) [-^) ' if']a:c--= ^ _ , 

[klJ := [l]a:c[2]a:c---[k]a:c. (15) 

The symbol {k}a:c should not be confused with < __2 r defined in (^), which 

looks similar but has a different meaning. The floor function [k\ denotes 
the largest integer < k. With these definitions the unnormalized A^-particle 
stationary state |A^) is built out of the vacuum by 

In addition, one has the left complete A^-particle state, defined as 

(iV| := (vac|^(5+)^. (17) 

■^Note that creation of particles corresponds to lowering a spin in the language of spin- 
operators. 
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For the left complete state it is straightforward to show that 

(iVl., =(iV-l|(l-nfc). (18) 

Prom exact solutions of systems with L < 4 a similar relation has been 
guessed for the unnormalized A^-particle state |A^): For all A: E {1, 2, . . . ,L} 
the following identity holds 

st\N) = {k}a-A^-nk)\N-l) (19) 

which can be proven by direct calculation. As a consequence one finds that 
the /-point density correlation functions obey the recursion relation 

{N\n,, . . . n.JiV) = {xi}a:c{N - l\n,, ... (1 - n,,)\N - 1). (20) 

With this recursion relation at hand it is easy to prove that 

P{ai,... ,afc-i,T,i,/?fc+2,-- - ,I3l) ^ /""^l (21) 

which suffices to show that the A^-particle states defined above are optimum 
ground states and hence stationary states. 

Note that B~ does not commute with H and so does not reflect a global 
symmetry. Both operators commute only, if they are restricted to the kernel 
of H. Therefore, the integrability found for a = c = g is most presumably 
lost if a / c. 



4. DENSITY PROFILES FROM RECURSION RELATIONS 

For the calculation of density profiles, in principle, the normalization 
of |A^) has to be computed. Unfortunately, a direct computation of this 
quantity is very involved and tedious. Therefore a step-by-step method has 
been chosen. Assume that Z{N — 1) is the normalization of — 1). Assume 
further that zn- is such that the normalization Z{N) of |A^) is given by 

Z{N) = znZ{N -I). (22) 

This gives with (pO|) 

Pk{N)= = {l-pk{N-l)). (23) 

Z (iV ) ZN 

In the A^-particle state naturally the identity Yl]^=i Pki^) = ^ holds. Thus, 

1 ^ 

ZN = j;^^{k}a:c{l-Pk{N-l)). (24) 

k=l 
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Fig. 1. Density profile of staggered hopping- model, a = 0.5, c = 2.5, L = 50 
and N = 10. The dashed line corresponds to analytical data from (p5|). The 
boxes represent Monte Carlo data. 



Finally, one finds a recursion relation for the local density of an A^- 
particle stationary state 

MN) = N fWl-/-^(^-l)) . (25) 

Ei=l{j}a:c(l-/>,(iV-l)) 

Starting from the normalized stationary state |vac), which belongs to = 
and has PkiO) = 0, the density profiles for all N can be computed. The results 
have been compared with Monte Carlo data (see Fig. |l| and Fig. ^) show- 
ing perfect agreement. Generically we find exponentially decaying density 
profiles. If one of the parameters is larger than one, whereas the other pa- 
rameter is smaller than one, the density profiles clearly reflect the sublattice 
structure of the Hamiltonian, since the slope of the profile is alternatingly 
positive and negative along the bonds (see Fig. |l|). The shape of the den- 
sity profile can be explained qualitatively in the following way: on the odd 
bonds hopping to the right is favoured because = 2 and a = 1/2. As a 
consequence, the density at the left end of an odd bond is smaller than the 
density at the right end. Along the even bonds one has the opposite scenario 
since c = 5/2 and = 2/5. Finally, the high density in the left part of the 
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Fig. 2. Density profile of staggered hopping-model. a = 2.0, c = 4.0, L = 50 
and = 10. The dashed line corresponds to analytical data from (p5|). The 
boxes represent Monte Carlo data. 

system is due to the fact that the maximal hopping amplitude can be found 
at the even bonds which favour hopping to the left. 

If both parameters a and c are larger or smaller than one, we find a 
different behaviour. Now the density profiles become monotonous functions 
of p and the sublattice structure is hidden (see Fig. |2[). 



5. A CAVEAT 

As presented above, optimum ground states come into play naturally 
if one encounters stochastic processes with reflective boundaries. Therefore 
one might conjecture that the stationary state of such processes is always of 
this special type. However, this is not true. In the following a Za-symmetric 
reaction-diffusion model will be presented, which has zero-energy ground 
states that are not (always) optimum ground states. 

The sites of a chain of length L are coloured with the elements of Z3, 
i.e., one has three types of particles (A; = 0, 1, 2) on a completely filled chain. 
Two particles on neighbouring sites react according to the following rules 
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(algebraic manipulations are done according to Z3) 

{k, k + 1) ^ {k + l,k) with rate a, 

(fc,fc + 2) ^ + + 1) with rate /3. 

Thus, the local Hamiltonian is given by 



(26) 
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(27) 



The operators act as ^ on sites j, j + 1 and as identity on the remaining 

sites of the chain, and hence the global stochastic Hamiltonian reads 



L-l 



(28) 



The exact solution for L = 2 shows a threefold degenerate ground state 

1^)2 = with A; e Z3. (29) 

These states are optimum ground states. 

Looking at L = 3 a new situation arises. One has again a threefold 
optimum ground state 



|V')3 = \k,k,k) with k G Z3. 
In addition, one finds two qualitatively different ground states 



(30) 



bi)3 = Yl 

^2)3 = Yl 



\k, k,k + 2) + ^\k,k + l,k + l) + -^\k, k + 2,k) 
a a + p 



\k,k + 2,k + 2) + ^\k,k,k + l) + -^\k, k + l,k) 
a a + p 



(31) 



These states are not optimum ground states. Moreover 

hi,2\^j)3 = -h2,d,Wj)d, / with j = 1,2. 



(32) 
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The occurrence of such states might be connected with the existence of an 
additional symmetry. 

The last result can also be understood in the context of MPA. As the 
two additional states are not optimum ground states, the cancelling vector 
X has to be nontrivial. Therefore, 

\<fj)3 = {W\V®V®V\V). (33) 

Using the fact that {W\ and \V) must be in the kernel of X one finds 

hi^2{W\V O P O V\V) = -{W\V A" V\V), 
h2,zW\V ®V® V\V) = {W\V (^X® V\V), 

and thus (|3^) holds. 

Most interestingly systems of length L = 4 and L = 5 - like the two- 
site system - possess only threefold degenerate optimum ground states. But 
then, for L = 6, two additional ground states show up again. This might 
be related to an additional symmetry occuring for system sizes which are an 
integer multiple of three. 

This example shows that reflective boundary conditions do not imply 
the existence of an OGS. However, it still leaves open the possibility that 
systems with reflecting boundary conditions and a unique stationary state 
have optimum ground states. 



6. CONCLUSIONS 

We have presented the exact solution for the stationary state of a 
stochastic process with staggered hopping rates. Using the similarity of 
the quantum formalism to the master equation the model could be inter- 
preted as a spin chain. Its ground state - corresponding to the stationary 
state of the stochastic process - was shown to be an optimum ground state, 
i.e. it is at the same time ground state of all local Hamiltonians which 
make up H. Furthermore the /-point density correlation functions have been 
determined recursively. 

Apart from the fact that the model is one of the very few exactly solv- 
able models with staggered interactions it might have practical a ppl ications. 
A similar model has been proposed by Kolomeiskv and Widomllll. It is a 
simplified ratchet model (for a recent review, see E3) for molecular r not ors 
effectively describing the motion of a motor protein on a micro-tubule 
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